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Abstract We report the results of a time-series CCD photometric survey of variable stars in 
the field of the open cluster NGC 2 1 26 . In about one square degree field covering the cluster, a 
total number of 21 variable candidates are detected during this survey, of which 16 are newly 
found. The periods, classifications and spectral types of 14 newly discovered variables are 
discussed, which consist of six eclipsing binaries systems, three pulsating variable stars, three 
long period variables, one RS CVn star, and one W UMa or S Scuti star In addition, there 
are two variable candidates, the properties of which cannot be determined in this paper. By 
a method based on fitting spectral energy distributions(SEDs) of stars with theoretical ones, 
the membership probabilities and the fundamental parameters of this cluster are determined. 
As a result, five variables are probably members of NGC 2126. The fundamental parameters 
of the this cluster are determined as: the metallicity to be 0.008 Zq, the age log(t)= 8.95, the 
distance modulus (to — M)o= 10.34 and the reddening value E{B — V)=Q.55 mag. 

Key words: open clusters: individual (NGC 2126)-stars: variables (general)-binaries: gen- 
eral 

1 mXRODUCTION 

Studies of variable stars in star clusters are crucial tools for understanding of stellar evolution and the nature 
of the host clusters. Recent works on variable stars in clusters can be found in Frandsen & Arentoft ( |1998l l; 
Park & Nemec ( 120001) ; Zloczewski et al. ( 120071 1 . Population I pulsating variables such as delta Scuti stars, 
7 Dor stars and slow pulsating B stars are populous in open clusters, and Population 11 pulsating variables 
such as RR Lyr-type stars and SX Phoenicis-type stars are rich in globular clusters. Stars in open clusters are 
not as crowed as globular clusters, so high-precision photometric results for open clusters can be obtained. 
Therefore, open clusters are very important targets to investigate pulsating variables. Eclipsing binary stars 
are also found in open clusters. For examples, 45 short-period eclipsing binaries were discovered in open 
cluster ColUnder 61 (Mazur et al I1995I I, seven eclipsing binaries are detected in the open cluster NGC 2099 
(Kang et al. I2007I I. sixty-one eclipsing systems were detected in the open cluster NGC 6791 (Marc hi et 
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The Beijing-Arizona-Taiwan-Connecticut (BATC) photometric system is a muhicolor sky survey. The 
system has shown the advantages for open cluster studies (Fan et al. 119961 Wu et al. 120051 120071) . This 
system has also contributes to search for variable stars in clusters. For instances, seventeen pulsating stars 
and ten eclipsing systems were found in open clusters NGC 188 (Zhang et al. l2002ll2004l) . Three pulsating 
variable stars were discovered in the direction toward NGC 4565 (Li et al. 120041) . Six eclipsing binaries 
and seven pulsating variables in the field of open cluster 2168 (Hu et al. 120051 ). In this paper, we search for 
variable stars in the field of NGC 2126. 

NGC 2126 (a2ooo=06''02™.55, ^2000= +49°52'; Z=163.23885, 6=13.13066° ) is a moderately rich typ- 
ical open cluster with several dozens of memberships scattered in a region of 5' ^ 6' in the constellation 
Auriga (Lynga[T987]). The first photometry was made in blue and red (AA 4300-6200 A) photographic ob- 
servation by Cuffey ( 119431 ), who estimated the distance to NGC 2126 as 950 pc. Sixty years later, the first 
CCD VRc Ic photometric observation was performed by Caspar et al. (120031 ). They found 6 variable stars 
in the cluster field, of which three could be members. They also obtained the fundamental parameters of 
the cluster: the age is 9.0<log(t)<9.4, the reddening value E{B - V) ^ 0'".2 ± O.^IS, and the distance 
d < 1.3 ± 0.6 kpc. To search for more cluster variable stars in the open cluster NGC 2126, we performed a 
dedicated CCD survey with the BATC system. In this paper, observations and data reduction are described 
in Section 2. Searching for variable stars is reported in Section 3. The physical parameters and the CMD of 
the Cluster are described in Section 4. Results and discussions are presented in Section 5. A brief summary 
is given in Section 6. 

2 OBSERVATIONS AND DATA REDUCTION 

2.1 Observations 

We observed NGC 2126 from 2004 December 4 to 2005 March 24. The observations were conducted with 
the 60/90 cm f/3 Schmidt telescope, which located at the Xinglong Station of the National Astronomical 
Observatories, Chinese Academy of Science (NAOC), with a 2Kx2K CCD camera at its main focus. The 
field view of the CCD covers about 1 deg^, with an image scale of 1.7" per pixel. The filter system consists 
of 15 intermediate-band filters with bandwidths 350 A-490 A, covering a spectral range 3000-lOOOOA, 
which are specifically designed to avoid most of the known bright and variable night-sky emission lines. 
In searching for variable stars observations, the/(5270A) and / (6660 A) band filters were employed. The 
exposure times were set as 300 seconds in /-band and 120 seconds in the /-band in each observation. In 
total, we got 487 frames of NGC 2126, 243 in/-band and 244 in /-band, respectively. 

2.2 Data Reduction 

Reduction of the CCD data including bias substraction and flat-fielding with dome flats was performed with 
an automatic procedure called PIPELINE I, which has been developed as a standard programme for BATC 
multicolor Sky Survey (Fan et al. 119961 ). The dome flat-field images were taken by using a diffuser plate in 
front of the correcting plate of the Schmidt telescope, the flat fielding technique has been verified (Wu et 
al. 120021 Zhou et al. 120041 ). The instrumental magnitudes of point sources in CCD frames were measured 
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Stetson ( 119871 ). The PIPELINE II program was performed on each single CCD frame to get magnitude 
of each point source by PSF (point spread function) fitting and the aperture photometries with different 
aperture. As a result, the list of PSF and aperture magnitudes in various filter bands were obtained for aU 
the objects detected. 

2.3 Flux Calibration 

In order to give the spectral types of the variables (see §3.1), we observed NGC 2126 with another 13 BATC 
filter bands (from c to p) in one night of semi-photometric conditions. Due to not being enough nights of 
photometric conditions, the flux calibrations of NGC 2126 are only defined in/-band and /-band by the four 
standard stars (Oke & Gunn[l983]l (BDh-174708, BDh-262606, HD19445 and HD849373). The flux of the 
four stars have been recalibrated by Fukugita et al. (1 19961 1. In this paper, the standard stars were observed 
between air masses of 1.0 and 2.0. The details of the flux calibration were described in Zhou et al. (I2001I I. 
Since the flux calibration of NGC 2126 in/-band and i-band have been defined by the four standard stars as 
described above, we can perform flux calibration of NGC 2126 in the other 1 1 bands based on the method 
called model calibration as presented by Zhou et al. (|1999l l. 

3 SEARCHING FOR VARIABLES 

In the frames, about 4416 stars are detected in the field of NGC 2126. To search for variable stars, we 
count on magnitude scatter of each measurement. The magnitude scatter of all stars in the /-band is shown 
in Figure 1. The abscissa is mean magnitude of each star. The dense parts are normal stars which are 
intrinsically stable in Figure 1. A star is selected as a variable candidate if 1) it presents a large deviation 
(> 3 rms) compared with that of the other stars of similar magnitude; 2) it has been measured on at least 
100 frames. Based on these conditions, hundreds of possible variable candidates were selected. Then, we 
checked the real-time light curves and the suspect time-series data of all candidates. The candidates showing 
noisy chaotic light variations, were rejected from the variable sample. After that, the number of the variable 
candidates decreased to about fifties. Furthermore, we checked the frames, and got rid of the stars that 
locate in the image edge, or saturated or are contaminated by the bright stars. Ultimately, a total of 19 
certain variable stars were identified from the candidates. Among these 19 stars, 5 variables named VI- 
V2,V4-V6 are known variables from Caspar et al. (I2003I I. the remaining 14 objects named V7-V20 being 
newly discovered. A known variable star named V3 in Caspar et al. ( I2003I I was not detected as a variable 
candidate in this paper We then check the data, and found that its variable amplitudes were too low to be 
detected. In addition, there are two candidates, the properties of which cannot be determined in this paper. 
The positions of all 21 variable candidiates in the field of NGC 2126 were marked in Fig. 2, and the newly 
found variables are marked in sequence of the distance away from the center of the cluster 

The analysis of the light-curves was made with the phase dispersion minimization (PDM) method 
(Stellingwerf |1978] l. In this method, the main characteristics of light curves such as the periods and am- 
plitudes of light variations, were given. 

According to the periods, spectral types and the behaviors of the light curves, the variable stars are 
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one eclipsing binaries with pulsating component, three long period variables, and three objects that cannot 
be confirmed based on the PDM method. The analyzed results, including their periods, variable amplitudes, 
mean magnitudes and spectral types, are also listed in Table 1. The real-time light curves of 6 variables are 
shown in Figure 3 . The phased light curves of V 1 , V2, V4, V5 and V6 are presented in Figure 4. The phased 
hght curves of 11 new variables are shown in Figure 5. The phased hght curves of V22 are shown in Figure 
6. 

4 THE PHYSICAL PARAMETERS AND CMD OF THE CLUSTER 

Based on the photometric data in the BATC 13 filter bands, the fundamental parameters of the open cluster 
NGC 2126 are determined by fitting the spectral energy distributions (SEDs) of stars in the cluster field 
with theoretical SEDs of Padova stellar evolutionary models (Girardi et al. 120001 l2002l Wu et al. l2005l l. At 
the same time, the membership probabilities of the stars including the variable stars studied in this paper, 
were determined based on the method of Wu et al. (120061 ). The best-fitting theoretical results of the cluster 
NGC 2126 are: the metalhcity Z=0.008Zq , the age log(t)=8.95, the distance modulus (m - M)o=10.34 
and the reddening value E{B — V) — 0.55 mag. The derived reddening is much large than that derived 
by Caspar et al. ( l2003t 

Malkov & Kilpio (12002) used various observational data to common interstellar extinction mod- 
els. They recommended to use the model proposed by Arennou et al. (11992) to calculate the extinction 
for object with distance great than 1 kpc and with galactic latitude \d\= 30°. Adopting the distance 
1.170 kpc of NGC 2126 calculated from above derived (m - M)o=10.34, the visual extincation Ay= 
1.467±0.369is determined based on the model of Arennou et al. ( 11992) . Adopting R^= 3.3, the extion- 
cation of NGC 2126 is estimated to be E{B - V)= 0.45±0.11, which is consistent with that derived in 
this study within the errors. 

Also based on the photometric data in the BATC 13 filter bands, the spectral types of the stars including 
the variable stars studied in this paper, can be given by fitting the BATC SEDs with the module SEDs of 
Gunn & Stryker (11983) . in which 74 spectral types from O to M6 and luminosities from main sequence to 
giant are included. 

The results of membership probabilities are shown in the last column in Table 1. In this table, "-" 
denotes that we did not obtain this parameter of this variable in this paper. The CMD of the observed field 
is shown in Figure 7. The filled circles denote the stars whose membership probabilities are more than 70 
percent. The black squares denote the variable stars. In addition, because VIO, V14, V16 and V19 do not lie 
e-band in the image, they are not plotted in the CMD diagram. The thick solid line denotes the theoretical 
isochrones with \og{t) = 8.95. The theoretical isochrones with the parameters derived by SEDs-fitting can 
fit star distributions in the CMD very well. The membership probabilities of VI and V2 are 0.91 and 0.65, 
respectively. They are on the main sequence of the cluster NGC 2126. Therefore, VI and V2 should be the 
members of NGC 2126. Considering the characteristics of variable stars and the positions in the CMD (Fig. 
7), V5 and V13 may be members of NGC 2126, although its membership probabilities are only 0.30 and 



Variable Stars in tlie Field of the Open Cluster NGC 2126 (I) 



5 



5 RESULTS AND DISCUSSIONS 

5.1 The known variables in NGC 2126 

There are a total of 5 known variable stars (Gaspar et al. I2003l lin the field of NGC 2126. Their basic 
parameters are listed in Table 1. During this observation, the good light curves for VI, V2, V4 ,V5 and V6 
were recorded. 

VI: Being similar to the previous work, we obtained two periods of VI, 0.82222 d and 1.64444 d 
comparing with 0.822235 d and 1.6447 d from Gaspar ( I2003I I. Gaspar ( I2003I I suggested that VI could be 
7 Dor or RS CVn star. The double-peak of the phase diagram is clear with the period of 1.64444 d (Fig. 
4). In our observation, we classified VI as A-type. However, the RS CVn criteria, which was defined by 
Hall |1976l and is widely used, suggests that the spectral type of RS CVn star is F^^K type. Kaye et al. ( |1999l l 
defined the 7 Dor phenomenon that consists of variable stars with an implied range in spectral type A7-F5 
and in luminosity class IV, IV- V or V. Therefore, VI may be 7 Dor star rather than a RS CVn. 

V2: It was discovered by Gaspar (I2003I I. and the periods around 0.5 d and 1 d were given, however 
neither period cannot give a continuous phase diagram with no gap. So, Gaspar ( 120031 ) did not determine 
the cause of the low-amplitude. In this paper, its spectral type is decided to be A-type. The phased light 
curves of Figure 4 show that V2 is becoming faint and its period may be longer than 2 day. 

V4: It is an Algol-type eclipsing binary. Gaspar (I2003I I observed only one minimum of this star and 
failed to determine its period. During our observations, we recorded three eclipse events. Using the PDM 
code, a probable period of 1.0966 d is derived. In addition, we determined its spectral type to be F-type. 

V5: Gaspar (I2003I I analyzed the frequencies of this star with Fourier analysis, and determined two 
frequencies /i =11.43 and /2=12.14 d^^ and /1//2 = 0.94. So they suggested V5 was a multiply 
periodic oscillations S Scuti star with non-radial modes. In this study, its spectral type is classified as A- 
type. Based on the PDM code, we derive a probable period to be 0.0827 d which is close to the /i, and the 
amplitude is 0.104 mag, very low. Combining the results of Gaspar (I2003I I and of this paper, V5 may be 
confirmed to be S Scuti star. 

V6: Gaspar (I2003I I showed that V6 was an eclipsing binary system with a S Set-type pulsating with 
Ppui — 0.129 d and Porb = 1-1732 d, and it might have fpui : forb = 1:9 resonance between the 
orbital motion and pulsation. Our results of Pp^i = 0.12957 d and Porb = 1.1732 d are consistent with 
Gaspar (I2003I I. Considering its spectral type of F-type determined in this paper, V6 may head Scuti star in 
an eclipsing binary. 

5.2 The newly discovered variables 

In this paper, sixteen new variables were discovered in the field of NGC 2126. According to the behaviors 
of the light curves, the spectral types and the period analysis, classifications were made. Among the sixteen 
new variables, four were classified as eclipsing systems, one as W UMa type star, one as EA star, one as W 
UMa star or 6 Scuti star, three as pulsating variable stars (one RR Lyr and two S Scuti), one RS CVn star, 
three long period variables, and two variable candidates that are not confirmed in this paper. Their basic 
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Fig. 3. The phased light curves of V7-V1 1, V14, V17-V21 are shown in Fig. 5. The phases light curves of 
V22 are plotted in Fig. 6. 

V7: Its position in the CMD , i.e. it is near the turn-off, shows, that it is a member of this cluster. 
It is difficult to determine its variable type. From the period of 0.395 d and the behaviors of the light 
curves, it looks like a W UMa binary system, which has a type of G or K and periods of 0.22 to 0.8 days 
(Binnendijk 1 1 965l l. However, according to the spectral type of A-type, and the Ught curves of sinusoidal 
variations, V7 may be a high amplitude 5 Scuti star The spectral types of 5 Scuti stars are from A to F with 
short periods of 0.02 to 0.3 days (Breger [T99T1) . 

V8: It is likely an eclipsing binary system with period of 0.28448 d. The phased light curves are plotted 
in Figure 5. 

V9: It posses typical W UMa-type light curves, and its spectral type is classified as a F-type, which 
consists with W UMa type, i.e. F^G. Therefore, V9 should be a W UMa eclipsing system with a period 
0.444462 d. 

VIO: This star may be a field eclipsing system, the period of which is 0.47743 d. The phased light curves 
are shown in Figure 5. 

Vll, V13 and V15: They may be field long period variable stars. As we know, the spectral types of red 
giants and supergiants are M, R, N or S, and the periods are up to 600 days. In our observations, Vll ,V13 
and V15 are all classified as M-types. Their real-time light curves also have long time variations (Fig. 3). 
Based on PDM code analysis, we obtained a probable period of Vll to be 38.5 d. The phase diagrams are 
plotted in Figure 5. The periods of V13 and V15 cannot be obtained in this paper. 

V14: It is a field Algol-type eclipsing system. Its hght curves have deep primary minima and flat bright 
parts. A period of 0.80042 d is obtained. 

V17: It is classified as A-type in this paper The light curves suggest it is likely a 5 Scuti variable with a 
period 0.1272 d. 

VI 8: About VI 8, its spectral type is classified as a K-type in this study. Analysis of the light curves 
reveals that it is likely an eclipsing binaries system with a period 0.20042 d. 

V19: Its real-time light curve posses the character of RR Lyra, asymmetrical, and increasing rapidly and 
decreasing slowly. So, V19 may be a RR Lyra system with P = 0.5375 d. 

V20: It is classified as an F - type star in this paper. From its light curves (Fig. 5), it may be a (5 Scuti 
star with P = 0.30315 d. 

V2 1 : The spectral type of V2 1 is classified as G-type in this paper. From §5 . 1 , we know that a RS C Vn 
is a later G ^ K type. So, V21 may be a RS CVn variable with a period 0.956 d. In / band, the light curves 
show that it is likely a RS CVn. However, in/band, they are much distorted. 

V22: The most confused object in this paper is V22. The analysis of light curves gives two periods, one 
is 0.22356 d and the other is 0.447124 d (Fig. 6). From the phased diagram of the period 0.22356 d, V22 is 
most likely an eclipsing contact binary system. On the other hand, the light curves of the period 0.447124 d 
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6 CONCLUSIONS 

In this paper, we presented a time-series CCD photometric survey of variable stars in the field of the open 
cluster NGC 2126. In this survey, 21 variables were discovered, including 5 previous ones (Caspar et 
al. l2003l l. Based on the SED fitting between the theoretical isochrones and the photometry, the membership 
probabilities of stars and the fundamental parameters of NGC 2126 were determined. The important results 
of this study are the following: 

1. We discovered 16 new variable stars in this survey: six were classified as eclipsing systems with four 
eclipsing binary systems (V8, VIO, V18 and V22), a W UMa type stai- (V9), an EA star (V14), three long 
period variables (VI 1, V13 and V15), three pulsating variable stars with a RR Lyr (V19) and two 5 Scuti 
(VI 7 and V20), and a RS CVn star (V21). About V7, its light curves show that it is most likely a W UMa 
type star. However, its spectral type is classified as A-type. So, it is very difficulty to confirm whether it is 
a W UMa star or 6 Scuti star. The more information will be needed in the future observation. In addition, 
two variable candidates(V12 and VI 6) were not confirmed in this paper. Among the 16 newly discovered 
variable stars, two (V7 and V13) may be the members of the open cluster NGC 2126. 

2. We also determined the mass of each variable star belonging to members of cluster NGC 2126, and 
the physical parameters of this cluster The masses of VI, V2, V5 and V13 are 1.1312, 0.9131, 1.4808 and 
O.4754Af0, respectively. The mass of VI suggests that it may locate at the blue edge of the instability strip; 
and its pulsations are driven by the modulation of the radiative flux from the interior of the star that is a 
relatively deep envelop convection zone (Guzik et al. 120001 Warner et al. l2003l) . As a 7 Dor, the variations 
of VI are consistent with the model of high order (n) and low spherical degree (I), nonradial, gravity (g)- 
mode oscillations (Kaye et al. ll999l l. Because we cannot determine the pulsating class of V2, its pulsating 
reason is unknown. Although the spectral type of V13 is known, the real-light curves are only a part of the 
whole, and the driving mechanism is not known. 

3. The physical parameters of the open cluster NGC 2126 are: the metallicity Z = O.OOSZq, the age 
log(i) = 8.95, the distance modulus (m - M)o = 10.34 and the E{B -V) = 0.55 mag. 
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Fig. 2 Observed CCD field (58' x 58') of the open cluster NGC 2126. Identification of variable 
stars are also marked. North is up and east is to the left 
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Fig. 3 Real-time light curves for six variable stars. 
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Fig. 4 The light curves for 5 known variable stars in /and i bands, respectively, with the periods 
marked in each frame. 
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Fig. 5 The light curves for new variable stars in / and i bands, respectively, with the periods 
marked in each frame. 




Fig. 6 The light curves for v22 with two possible periods in i bands. 



Variable Stars in tlie Field of the Open Cluster NGC 2126 (I) 




Fig. 7 The color-magnitude diagrams of NGC 2126. The solid line show the log(i)=8.96 theorit- 
ical isochrone (Girardi et al. 120001) . The filled circles denote stars observed whose probabilities 
to be members are more than 70 percent. The black squares denote the variabel stars. 
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Table 1 Observation summary 21 variable star in open cluster NGC 2126. 



ID 


• 

X cnoQ 


012000 


t'2000 


_ 

TTlil ) 


error 




opec 


VallaUlC Ly^c 


Mem 


VI 


1.64444 


06:01:44.09 


+49:56:30.3 


13.868 


0.02 


0.198 


A 


7 Dor 


0.91 




or 0.82222 


















V2 




06:01:57.36 


+49:58:54.8 


14.356 


0.025 


0.117 


A 


" 


0.65 


V4 


1.0966 


06:02:21.28 


+49:52:37.0 


15.507 


0.091 


0.637 


F 


EA 


0.00 


V5 


0.0827 


06:02:26.38 


+49:51:56.5 


12.930 


0.05 


0.104 


A 


<5Sct 


0.30 


V6 


1.1732 


06:02:37.99 


+49:53:02.5 


14.211 


0.045 


0.479 


F 


EA+<5Sct 


0.00 




or 0.129572 


















V7 


0.395 


06:04:04.14 


+49:45:53.4 


12.671 


0.066 


0.240 


A 


W UMa or 5 Scuti 




V8 


0.28448 


06:01:14.33 


+49:40:58.9 


16.364 


0.184 


0.808 




eclipsing binaries 


0.00 


V9 


0.444462 


06:03:34.89 


+49:36:45.0 


15.626 


0.083 


0.454 


F 


WUMa 


0.00 


VIO 


0.47743 


06:00:46.04 


+50:00:34.9 


17.032 


0.163 


0.728 




eclipsing binaries 


0.00 


Vll 


38.5 


06:04:34.69 


+49:58:33.4 


11.032 


0.045 


0.170 


M 


LPV 


0.00 


V12 




06:01:55.94 


+49:32:08.5 


10.122 


0.024 


0.134 


- 


uncertain 


0.00 


V13 


- 


06:02:21.61 


+50:02:35.5 


11.007 


0.018 


0.090 


M 


LPV 


0.46 


V14 


0.80042 


06:00:22.78 


+50:03:40.4 


14.013 


0.058 


0.392 




EA 




V15 




06:03:07.49 


+50:18:13.8 


11.395 


0.017 


0.111 


M 


LPV 


0.00 


V16 




06:00:32.51 


+49:34:27.4 


11.692 


0.012 


0.085 




uncertain 




V17 


0.0636 


06:04:26.43 


+49:27:39.5 


11.494 


0.012 


0.093 


A 


S Scuti 


0.00 


V18 


0.2004 


06:05:46.45 


+49:42:11.0 


16.084 


0.160 


0.721 


K 


eclipsing binaries 




V19 


0.5375 


06:00:07.54 


+49:30:07.5 


14.640 


0.042 


0.279 




RRaLyr 




V20 


0.30215 


06:04:39.65 


+50:20:25.8 


15.228 


0.111 


0.709 


F 


5 Scuti 




V21 


0.956 


06:05:50.90 


+50:18:01.8 


12.961 


0.022 


0.169 


G 


RS CVn 


0.02 


V22 


0.22356 


06:05:42.28 


+50:20:27.0 


16.642 


0.201 


1.158 




eclipsing system 


0.00 



or 0.447124 




C\2 



a; 
m 

Oh 



d 



9'tT 




C\2 



9T S'9T LI Q'^T 



a; 
m 

Oh 



d 



81 




C\2 



91 9'9T L\ Q'^T 



81 



a; 



d 



CD ?H 



-a 

CO 
CO 

d 



03 



— •••• 



C\2 



• • • * 



4r. 



: • 

•••••• 

%• •> • 
••• # 

• ••• 

• • • 

J I I I •! 



91 9'9T 



a; 



d 





1 1 1 1 1 1 1 
• 


1 1 1 1 1 1 1 














m 








^ I 


i 












• 1 — 1 




♦ 




_+ 


























• 








• 






in 


• 






to 


• • 


+ ¥ 




CO 


• 


+ 






m 








i 






> 


1 


♦ 






















• 








• 








• 


i. 






1 • * , , 1 , 


1 1 1 1 1 1 -t-hnri 





o 



IT 9"TT Z\ Q'2T 




I I \ I \ I \ I \ I \ U 



